Optimal Control of Stochastic Magnetization Dynamics by Spin Current 
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Fluctuation-induced stochastic magnetization dynamics plays an important role in magnetic 
recording and writing. Here we propose that the magnetization dynamics can be optimally con- 
trolled by the spin current to minimize or maximize the Freidlin-Wentzell action functional of the 
system hence to increase or decrease the happening probability of the rare event. We apply this 
method to study thermal-driven magnetization switching problem and to demonstrate the advantage 
of the optimal control strategy. 



One of the key issues in the modern magnetoelectron- 
ics is the reliable control of the magnetization dynam- 
ics. Compared with the conventional Oersted field gen- 
erated by the electric current, the spin transfer torque 
(STT) [H, 0| acting on the magnet by the spin current has 
been proven to be a more efficient manipulation method 
as the magnetic structures are miniaturized towards the 
nanoscalc S, 4,)- Various devices based on this principle 
has been realized, such as the STT-magnetoresistive ran- 
dom access memory, spin-torque oscillators, spin logic 
etc[i|. Although the effect of STT on the magnetization 
dynamics has been extensively studied, its effect on the 
magnetization fluctuation, which can have a prominent 
impact on the performance of the spintronics devices 0- 
[lol |. has not received enough attentions. Recently, it was 
demonstrated experimentally that spin current can either 



suppress or enhance the magnetization fluctuations [11 
which coincides with the theoretical studies (l2. 13 1. Fur- 
thermore, it raises the interest to control the magneti- 
zation fluctuation by spin current and different control 
strategies have been proposed[l3, 15|. However, there 
still exists the controversy whether the effective damping 
should be decreased or increased in order to suppress the 
magnetization fluctuation IJ, [15[ , and the choice of con- 
trol strategies is somewhat arbitrary without solid the- 
oretical foundation . Thus deeper understanding the 
effect of spin current on the stochastic magnetization dy- 
namics and searching for better control strategy are de- 
manded for further applications of the spin current in 
spintronics devices. In this letter, we will analyze how 
the spin current can affect the stochastic magnetization 
dynamics based on the large deviation principle, and how 
to optimize the spin current pulses to control the mag- 
netization fluctuation. The application of the idea to 
the problem of thermal-driven magnetization switching 
in the presence of STT will be demonstrated as an ex- 
ample. 

The thermal-driven magnetization dynamics under 
STT is usually given by the stochastic LLG equation|12 
in the case that the quantum effect is negligible For 
given initial magnetization configuration, the fluctuating 
magnetic field due to the thermal noise will make the 
trajectories of the magnetzation configuration randomly 
deviate from its deterministic trajectory when there is 



no thermal noise. When the noise amplitude is small, 
most of these random trajectories are around the deter- 
ministic trajectory, and only few trajectories can deviate 
from the deterministic trajectory largely. Thus the mag- 
netization dynamics at long-time scale can be regarded 
as the combination of the quasi-deterministic motion for 
most of the time and some large deviation trajectories 
happened rarely and randomly. Nevertheless, it is be- 
cause the stochastic magnetization dynamics is mainly 
governed by these rare trajectories that the control of 
the stochastic dynamics is reduced to the control of these 
rare trajectories. 

For simple, we will consider the dynamics of the single- 
domain magnet with magnetization vector M = Msin 
and volume V, where Ms denotes the constant magneti- 
zation magnitude and m denotes the unit direction vec- 
tor. The case for multidomain magnetic structures is 
more complicated but can still be treated in the simi- 
lar way. Generally, the happening probabilities P[m(i)] 
of the rare magnetization trajectories {m{t)} in stochas- 
tic magnetization dynamics satisfy the large deviation 
principle Til- 19j| . i.e. taking the exponential form 



P[m(i)] 



,-S[m(t),7=(t)]/e 



(1) 



where the Freidlin-Wentzell (FW) action functional 
S[in{t),I^{t)] depends on the applied spin current 
and the parameter e reflects the noise amplitude. In or- 
der to increase (or decrease) the happening probability of 
certain magnetization trajectory {m(t)}, one can adjust 
the applied spin current I" (t) to decrease (or increase) the 
FW action functional. Specially, the maximal (or min- 
imal) happening probability for {m(t)} is achieved by 
minimizing (or maximizing) the FW action functional. 
Usually, certain constrait conditions for the spin current 
should be set up, such as the total consumed energy dur- 
ing the control should be constant, etc. Then the optimal 
spin current I^ptit) to control the happening probability 
of the given magnetization trajectory {m(i)} is given by 
the variational problem f20t 



5{S[xn{t),P{1^]+\F[Pm^Q, 



(2) 



where we have set the constraint conditions F[I'^{t)] — 
for the spin current, and A is the Lagrange multiplier. 
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The dynamics of the single-domain magnet are driven 
by the deterministic effective magnetic field given by the 
micromagnetic energy density, the STT due to the spin 
current, and the stochastic magnetic field arising from 
the thermal fluctuation. Explicitly, the effective mag- 
netic field is He// = — V „!-£'( m) /Ms, with the micro- 
magnetic energy density £'(m). The spin current is 
defined by its amplitude aj and spin-polarization vector 
P, and we will only consider the adiabatic STt[i|, 13] al- 
thou gh t he non-adiabatic STT can be important in some 
cases I21j. The fiuctuating magnetic field is assume as 
-^eW, where the amplitude e = ^"j^^^ is proportional to 

temperature T and Gilbert damping coefficient a, and W 
is the Gaussian white noise process satisfying (W{t)) = 
and {Wi{t)Wj{t')) = 5ij5{t—t'). Here we have introduced 
the Boltzmann constant fc^ and the gyromagnetic ratio 
7. Then the stochastic LLG equation can be written in 
the compact formfiot 

m = b(m) + ^/icr(m)W, (3) 

where the vector b and the matrix a are defined as 



b(m) 
cr(m) 

Here we have 7' — 

symmetric matrix Ka and symmetric matrix Ks are[lt! 



^ 7, and the elements of the anti- 



l+a 



(Ka. 



{Ks) 



In this case. 



the expression of the FW action functional S[ra.{t) , P [t)] 
will be[il [H 



(4) 



where ti and tf are the initial and final time for the tra- 
jectory {m}. Equations © and (0]) are the theoretical 
foundations for us to optimally control the targeted tra- 
jectory {m} of the sinlge-domain magnet with spin cur- 
rent I" . 

For the above optimal control problem, the usual situ- 
ation is that many targeted trajectories exist. For exam- 
ple, the magnet may switch from one stable configuration 
to another via different paths. Among these targeted tra- 
jectories, the one which minimizes the FW action func- 
tional S'[m(i), J*] for fixed spin current P would be the 
most probable trajectory, which is the so-called optimal 
pa</i[20|. The optimal path provides the best prediction 
for the random and rare trajectories happened in prac- 
tice, and the corresponding happening probability is P ^ 
g-y[m„pt,/1/e^ where V\mopt,P] = min{S'[m(i), P]} is 
the quasi-potential[l^ Thus, the optimal control 

problem is further reduced to find the optimal spin cur- 
rent to minimize or maximize the quasi-potential in 
order to increase or decrease the happening probability 
of the optimal path. 

In principle, the double optimization problem formu- 
lated above can be solved by various optimization al- 
gorithms. One special but important case is that the 



spin current P is weak and can be treated perturba- 
tively in the magnetization dynamics. Then the optimal 
path {niopt} can be assumed to be independent on the 
spin current, and the quasi-potential V[raoptiP] can be 
approximated to the first order of aj as ^[mopt,/''] = 
Vo[mopt] + Ay[mopt, P], where 



|a-i(m. 



opt 



(5) 



= ^ ' aj7'(a-i(m„pt - bo), a-'KsP)dt. (6) 

For convenience, we have introduced the vector bo = 
aHeff, and notice that a, Ks, He// are all functions of 
niopt- Thus, the presence of a weak spin current P gives 
the change of the initial quasi-potential Vq in the amount 
of A.V, and the happening probability of the optimal path 
{m.opt{t)} has been exponentially changed by the factor 
g-AV/e^ The optimal spin current /q^j is obtained by 
minimizing or maximizing AV , i.e. 



S{AV[moptit), P{t)] + AP[P(t)]} = 0. 



(7) 



Thermal-driven magnetization switching is a typical 
example of the stochastic dynamics described above, and 
the concept of "effective temperature" has been intro- 
duced to the Neel-Brown law to take account into the 
effect of spin current on the magnetization switching 
probability?-'^ This fact can be easily verified from 
Eq. ([5]) and ([B]) for a weak spin current with constant 
amplitude, where Ay is now proportional to aj and the 
magnetization switching probability is formally written 
as P ^ g~(i~^)^o/'^_ Here, Qc is the critical value for the 
spin current to switch the magnet, and is given by 



a- 



Jl^\a-^KsIl,fffdt 
'Jl'{a-^KsII,ff,a~^KsP)dt 



(8) 



To get Eq. ([S]), we have used the fact that in the ab- 
sence of STT, the optimal path {niopj} satisfies the 
eauation[l8| 



eff- 



(9) 



From Eq. ([5]) , it is seen that the value of Uc is dependent 
on the choice of the spin polarization vector P of the 
spin current. Especially, if P is always in anti-parallel 
with the effective magnetic field Hg//, then Oc will be 
negative and the magnetization switching probability P 
will be decreased exponentially. From Eq. (j3|) , one might 
regard that the Gilbert damping should be effectively 
enhanced by the spin current in order to stabilize the 
magnet |_15j. One the other hand, P should be parallel 
with He// in order to increase P. Furthermore, Eq. (|8]) 
gives a rough estimation of the critical spin current mag- 
nitude as flc ~ amax{|He//|}, which is proportional to 
the Gilbert damping coefficient and the maximal effective 
fieldfU. 
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If the spin current amplitude is assumed to be time- 
dependent, the optimal pulse shape a'j'^{t) to control the 
optimal path {rriopt} can be obtained from Eq. ([7]). One 
natural constraint condition would be that the total dis- 
sipation energy of the spin current pulse should be con- 
stant, i.e. ^ a^j{t)dt = £. Then the general solution 



opt 



(t) of Eq. is given by 



opt 



it) 



A 



opt/ 



0. (10) 



and the Lagrange multiplier A is determined from the 
constraint condition as 



-'KsB.,ff,a-'KsPoptfdt]^. (11) 



This spin current pulse gives the change of the quasi- 
potential as AV — —2\£. According to Eq. pO|) . the 
sign of A is postivie (or negative) when the optimal spin 
polarization vector Popt is parallel (or antiparallel) with 
the effective magnetic field Hef / , which will enhance (or 
suppress) the happening probability of the optimal path 
mopt(i) by the factor e^^^l"^ . 

Based on the theoritical analysis above, we use the 
stochastic LLG equation to simulate the optimal control 
of the thermal-driven magnetization switching by spin 
current numerically. Considering a ferromagnetism film 
with easy axis along z-axis, and demagnetization field 
direction along x-axis, then the energy density E(va) is 



given as 



(12) 



Here, we assume that no external magnetic field is ap- 
plied, the anisotropic field is Ha = 0.05 T, and the de- 
magnetization field is AttMs = 1.2 T. Besides, we set the 
Gilbert damping coefficient a = 0.03, the temperature 
T = 300 K, and the magnet volume V = 1500 nm^. The 
initial direction of the magnet is assumed as m = (0, 0, 1), 
and the stochastic LLG equation ^ is simulated with 
the 4th-order Runge-Kutta method, where the time step 
is set as 1 ps. We also got the amplitude of the critical 
spin current numerically as Oc = 0.024 at zero tempera- 
ture. Notice that we have chosen a relatively small mag- 
net otherwise the switching time will be too long to be 
numerically simulated. 

The energy profile £'(m) has two saddle points, i.e. 
(0, 1, 0) and (0, —1, 0), thus there exist two corresponding 
optimal paths for the magnet to switch from the initial 
stable point (0,0,1) to the other (0,0,-1). Fig. ^a) 
shows the optimal path {ra^^^{t)} from (0,0,1) to 
(0, 1,0), which is obained from Eq. As one can see, 
there are some oscillations around the stable point before 
the magnet arrives at the saddle point, and this process 
takes about 2 ns. For comparision. Fig. [Ijb) gives one 
random switching trajectory of the magnet simulated by 
the stochastic LLG equation ([3]) with no spin current. 



It is found that the main feature of the real switching 
trajectories are indeed qualitatively captured by the op- 
timal path, although there are some quanlitative differ- 
ences because the noise amplitude is not small enough 
in our simulations. We can also see that the magnet 
randomly moves around the stable point for most of the 
time (about 1249 ns in this sample trajectory), and the 
switching event happens rather rarely. 

The optimal spin current pulse I^^{t) to enhance 



the happening probability of the optimal switching path 
{m^pj(t)} are given in Fig. [TJc) and (d). The absolute 
value of a°f*{t) is determined by the constraint condi- 
tion, and its shape shows that the best control strategy 
is to allocate the spin current when the magnet has de- 
viated from the stable point significantly. The optimal 
spin polarization direction is set to be in parallel 
with the effective magnetic field Heff- In contrast, the 
optimal spin current pulse to suppress the optimal path 
{m^pj(t)} should have an opposite spin polarization di- 
rection. For the energy profile considered here. He/ f 
has no y-component, and its a;-component can be quite 
large even for a small m°P* due to the large demagneti- 
zation field. 
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FIG. 1. (Color online) (a) optimal path {m^pj(t)} for the 
magnet to switch from (0, 0, 1) to (0, 1, 0); (b) a sample mag- 
netization switching trajectory simulated by the stochastic 
LLG equation; (c) the amplitude of the optimal spin current 
pulse corresponding to {in^p((t)}; (d) the spin polariza- 
tion vector of the optimal spin current pulse to enhance 
the happening probability of {m^p^{t)}. 

As a test, we calculated the mean switching time 
(MST) of the magnet under different control strategies, 
where 10^ random trajectories are generated for each 
case. First, the MST without spin current is about 
567 ns. Then a constant spin current with aj = 0.05ac 
will increase the MST to be about 913 ns if P = 
(0,0,-1), and decrease the MST to be about 353 ns if 
P = (0,0,1). While for the optimal control, the cor- 
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responding spin current pulse 1^^^ is applied only if the 
magnet randomly arrives at the critical position m that 
= nic- We have chosen — 0.89,0.81,0.65 seper- 
ately, and the initial point of the corresponding spin cur- 
rent pulses are denoted as 1, 2, 3 in Fig. H^c). The ampli- 
tudes a°f* are given such that each pulse will contain the 
same energy as a 2-ns constant pulse with aj = O.OSoc. 
Thus, the number of pulses used directly reflects the en- 
ergy consumption during the control process. 



TABLE I. Mean switching times r (unit : ns) and number of 
spin current pulses n when the switching probability is either 
suppressed (s) or enhanced (e). rric : the critical values to 
activate the control spin current pulse. Subscript a : the 
optimal spin current pulses I opt is applied; subscript b : the 
spin polarization vector is fixed as (0,0,-1) or (0,0,1) for 
suppressing or enhancing the switching respectively. 











0.89 


517/430 1067/892 


28/26 


292/249 


0.81 


739/263 1073/390 


32/13 


297/110 


0.65 


515/32 903/57 


67/4 


347/22 



The resulting MSTs r and mean pulse numbers n are 
shown in Table HI for both the cases to suppress and en- 
hance the switching probability. Unfortunately, we found 
that the optimal pulses didn't increase the MSTs when 
rric = 0.89 and 0.65. For rric — 0.81, the situation is bet- 
ter since the MST has been increased to 739 ns with 
263 pulses, but it still has no obvious advantage com- 
pared with the constant control current case. The reason 
of the failure is that the noise amplitude is not small 
enough and the real switching paths m(t) can deviate 
from the optimal path significantly so that the control 
pulse will in fact increase their happening probabilities. 
To eliminate this effect, we have set the spin-polarization 
direction as P = (0, 0,-1), and the results indeed become 
much better, as shown by and in Table H Especially 
when rric = 0.65, the MST is 903 ns which is close to 
the constant spin current case, but the number of pulses 
nl has been greatly decreased to 57. This shows that 
the optimal control theory is helpful to find the way to 
save energy even when the noise amplitude is not small 
enough. 

The advantage of the optimal control strategy is fully 
verified for the case to increase the switching probability. 
In TablelH we found that the MST is decreased to 28 ns 
with only 26 pulses for = 0.89. If = 0.65, the MST 

= 66 ns is larger, but merely 4 pulses in average is used 
here. For comprison, we also list the results and 
if P is set as (0, 0, 1) for the spin current pulses, which 
are only a little better than the constant spin current 
case. Obviously, the optimal control strategy successfuly 
increased the switching probability with greatly reduced 
energy consumption. 

In conclusion, we have discussed how to use the spin 



current to optimally control the stochastic magnetiza- 
tion dynamics based on the large deviation principle, and 
shown the advantage of the proposed control strategy by 
applying it to the thermal-driven magnetization switch- 
ing problem. The happening probabilities of the rare 
events, which are the dominant part of the stochastic 
magnetization dynamics, are controlled by changing the 
corresponding quasi-potential with optimized spin cur- 
rent pulse. Further generalization of the idea to the cases 
including the multi-domain magnetization configuration, 
quantum fluctuations etc. are necessary for more practi- 
cal applications . 



Y. W. thanks Dr. G. D. Chaves-O'Flynn for helpful 
discussions. This work is supported by the Hong Kong 
University Grant Council (AoE/P-04/08). 



[1] J. C. Slonczewski, J. Magn. Magn. Mater. 159, LI 
(1996). 

L. Berger, Phys. Rev. B 54, 9353 (1996). 
D. C. Ralph and M. D. Stiles, J. Magn. Magn. Mater. 
320, 1190 (2008). 

A. Brataas, A.D. Kent, and H. Ohno, Nature Materials 
11, 372 (2012). 

J. W. Lau and J. M. Shaw, J. Phys. D: Appl. Phys. 44, 
303001 (2011). 

D. Waller and A. Moser, IEEE Trans. Magn. 35, 4423 
(1999). 

S. Urazhdin, N. O. Birge, W. P. Pratt, Jr., and J. Bass, 
Phys. Rev. Lett. 91, 146803 (2003). 

R. H. Koch, J. A. Katine, and J. Z. Sun, Phys. Rev. Lett. 
92, 088302 (2004). 

I. N. Krivorotov, N. C. Emley, A.G. F. Garcia, J.C. 
Sankey, S. I. Kiselev, D. C. Ralph, and R. A. Buhrman, 
Phys. Rev. Lett. 93, 166603 (2004). 

X. Cheng, C. T. Boone, J. Zhu, and I. N. Krivorotov, 
Phys. Rev. Lett. 105, 047202 (2010). 
V. E. Demidov, S. Urazhdin, E. R. J. Edwards, M. D. 
Stiles, R. D. McMichael, and S. O. Demokritov, Phys. 
Rev. Lett. 107, 107204 (2011). 

Z. Li and S. Zhang, Phys. Rev. B 69, 134416 (2004). 
Y. Wang and L. J. Sham, Phys. Rev. B 85, 092403 
(2012). 

M. Covington, U. S. Patent No. 7,042,685 (9 May 2006). 
S. Bandopadhyay, A. Brataas, and G.E.W. Bauer, Appl. 
Phys. Lett. 98, 083110 (2011). 

M. Freidlin and A. Wentzell, Random Perturbations of 
Dynamical Systems (Springer- Verlag, New York, 1998). 
H. Touchette, Phys. Rep. 478, 1 (2009). 
R. V. Kohn, M. G. Reznikolf, and E. Vanden-Eijnden, J. 
Nonlinear Sci. 15, 223 (2005). 

G. D. Chaves-O'Flynn, D. L. Stein, A. D. Kent, and E. 
Vanden-Eijnden, J. Appl. Phys. 109, 07C918 (2011). 
V. N. Smelyanskiy and M. L Dykman, Phys. Rev. E. 55, 
2516 (1997). 

S. Zhang and Z. Li, Phys. Rev. Lett. 93, 127204 (2004). 



[2 

[3: 
[4: 
[5: 
[6: 
[t: 
[s: 
[9: 

[10 

[11 

[12 
[13 

[14 
[15 

[16 

[17 
[18 

[19 

[20 

[21 



